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ABSTRACT 


The  volatile  radioactive  species  produced  by  irradiation  of 
AON  -  201  Reactor  fuel  elements  were  obtained  by  irradiation  of  the  fuel 
in  an  evacuated  inert  container.   The  following  11  constituent  nuclides, 

including  two  radioactive  daughters,  were  identified  by  gamma  ray  spee- 

.131  v  131  v  131m  v  133  v  135  v  138  „  85m  v   8?  v   88 
trometry:   I   ,  Xe   ,  Xe    ,  Xe   ,  Xe   ,  Xe   ,  Kr   ,  Kr   ,  Kr   , 

I  OQ  QQ 

Cs   ,  and  Rb   .   Additionally,  the  possible  presence  of  the  following 

.  , .       T133  T13A  T135   r136  v  137  „  88  „  89    . 
species  was  indicated:   I   ,1   ,1   ,1   ,  Xe   ,  Kr   ,  K   ,  and 

89 
Rb   .   The  spectra  studied  were  taken  with  a  Nal  (Th)  crystal  and  they 

contained  34  photopeaks.   An  automated  computerized  method  for  half  life 

analysis  of  complex  photopeaks  is  described. 
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ANALYSIS  OF  VOLATILE  RFACTOR  FUEL  DECOMPOSITION 
PRODUCTS  BY  SPECTRAL  ANALYSIS 

I.   Introduction 

The  Naval  Postgraduate  School  Reactor  Facility  contains  an  AON-201 
nuclear  reactor  that  has  been  modified  to  produce  up  to  1000  watts  for 
short  periods  of  operation.   The  core  and  rods  of  the  reactor  contain 
uranium  dioxide  suspended  uniformly  in  a  polyethylene  matrix  (1) .   During 
operation,  of  course,  numerous  fission  products  are  formed,  and  the 
polyethylene  is  radiolyticallv  decomposed  by  recoil  and  the  conglomerate 
photon  and  particulate  radiation  present.   Many  of  the  radioactive  sub- 
stances thus  formed  are  volatile  so  they  must  be  withdrawn  from  the 
reactor  in  a  controlled  manner,  and  this  is  done  routinely.   Since  the 
scheme  of  the  disposal  is  to  exhaust  these  substances  into  the  atmosphere, 
the  characterization  of  these  radioactive  substances  has  considerable 
practical  application.   Also,  the  suspension  of  the  uranium  dioxide  in 
polyethylene  in  the  reactor  is  a  system  which  has,  potentially,  a  very 
interesting  chemistry,  because  of  the  multitude  of  organic  and  inorganic 
hot  atoms  that  are  mixed  intimatelv  together. 

Much  is  known  about  fission  products,  per  se,  and  radiation  effects 
on  polyethylene,  and  some  research  has  been  done  on  the  gaseous  effluents 
at  the  NPS  reactor.   The  composition  of  the  gaseous  effluents  from  the 
NPS  reactor  has  been  reported  on  by  Orissom  (2),  Bredderman  (3),  and 
Laskaris  (4) .   The  first  two  investigators  assumed  that  almost  all  the 
volume  of  gas  generated  was  hydrogen.   They  both  calculated  the  specific 
activity  of  the  core  gas  and  Bredderman  identified  six  radioactive  nuclides 
by  gamma  ray  spectrometry.   Bredderman  used  the  capable,  but  laborious, 
method  of  Heath  (5)  to  determine  photopeak  half  lives.   Laskaris  used 
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simultaneous  gas  chromatography  and  scintillation  counting  to  show  that 
reactor  gas  was  approximately  50%  hydrogen  by  volume,  and  that  xenon  and 
krypton  were  present. 

The  previous  in  house  efforts  on  the  analysis  of  reactor  effluents 
made  considerable  progress  toward  explaining  the  gross  characteristics  of 
the  situation,  however,  they  were  handicapped  in  several  ways.   Their 
samples  were  taken  at  the  exhaust  mainfold  of  the  reactor,  and  this  is 
separated  from  the  reactor  core  by  about  4  meters  of  5  mm  I.D.  copper 
tubing.   It  is  significant  in  this  regard  that  the  cylindrical  core  con- 
tains, in  addition  to  the  fuel,  an  annular  ring  of  graphite  that  is  about 
2.5  centimeters  thick,  25  centimeters  in  diameter,  and  25  centimeters 
high.   The  fuel  -  moderator  system  referred  to  earlier  consists  of  about 
3  kilograms  of  uranium  dioxide  and  about  12  kilograms  of  polyethylene. 
Thus,  samples  taken  at  the  manifold  may  not  have  been  truly  representative 
of  the  gases  produced  in  the  core.   Reactive  species  present  in  only  small 
quantities,  for  example,  iodides  and  bromides,  have  a  very  high  proba- 
bility of  being  totally  absorbed  before  reaching  the  sample  bottle.   Further- 
more, every  sample  was  almost  certainly  different,  because  it  contained  an 
unknown  amount  of  aged  gas.   The  composition  of  this  aged  gas  was  depen- 
dent on  both  the  length  of  time  the  active  constituents  had  decayed,  and 
the  relative  rates  of  diffusion  of  all  the  constituents  out  of  the  poly- 
ethylene mass.   Another  uncertainty  was  introduced  by  the  fact  that  the 
core  is  pressurized  with  water  pumped  nitrogen  which  contains  greater  than 
atmospheric  concentrations  of  most  of  the  rare  gases. 

It  was  believed  that  the  foregoing  handicaps  were  sufficiently 
important  that,  if  they  could  be  overcome,  another  attempt  at  analysis  of 
the  core  gas  was  warranted.   The  problems  were  circumvented  and  the  analysis 
was  Started •   The  vital  factor  In  meeting  the  problems  was  the  decision 


to  obtain  fuei  gas  in  a  much  more  controlled  manner  by  irradiation  or 
reactor  fuel  in  an  evacuated  inert  container.   This  technique  changes  the 
project  from  analysis  of  the  aforementioned  core  gas  to  analysis  of  fuel 
gas.   The  principal  implications  of  this  are  that  the  samples  now  analyzed 
will  be  void  initially  of  the  long  lived  constituents  of  the  aged  gas, 
but  it  will  be  possible  to  detect  many  of  their  short  lived  precursors. 
The  most  important  benefits  to  be  derived,  however,  are  that  the  samples 
will  be  reproducible  and  will  contain  greater  auantities  of  reactive  and 
trace  constituents. 

Additionally,  the  project  would  afford  an  excellent  opportunity  to 
explore  methods  of  applying  the  two  recently  acquired  single  channel 
analyzers,  with  rapid  digital  read  out,  to  the  analysis  of  complex  gamma 
spectra.  It  was  felt  that  there  was  considerable  opportunity  to  contri- 
bute to  the  state  of  the  art  of  unraveling  complex  gamma  spectra  by  auto- 
mating and  computerizing  data  acquisition  and  reduction.  This  was  taken 
as  a  secondary  objective  of  this  research. 
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SAMPLING  PROCEDURE:   Evacuate  entire  system  to  less  than  1  micron  of  Hg; 
Seal  at  F  and  remove  G;  Irradiate:  Pump  gas  into  ampul  with  liquid 
nitrogen;  Seal  at  D  and  remove  ampul  for  counting. 

COMPONENTS:   A,  Al  container;  B,  SS  bellows  valve;  C,  SS  ST;  E,  quartz 
tubing;  G,  S J . 


Figure  1 

II.   Experimental 

Reactor  fuel  elements  were  irradiated  in  a  hermetically  sealed 
iluminum  container  in  an  access  part  of  the  reactor,  and  a  sample  of  the 
gas  pumped  with  liquid  nitrogen  directly  to  a  sealed  quartz  amput  for 
gamma  ray  spectrum  analysis.   Quartz  was  used  to  preclude  spurious  peaks 
from  activation  of  the  ampul.   An  attempt  to  activate  a  sample  of  quartz 
stock  showed  it  produced  no  interfering  peaks.   Normally  one  fuel  ele- 
ment was  irradiated  to  produce  the  gas  sample,  and  a  typical  irradiation 

9     2 
was  10  n/cm  /sec  for  15  minutes.   A  fuel  element  is  a  circular  disc 


4  cm  in  diameter  and  1  cm  thick  that  contains  1.0  gram  of  U 


235 


The 


uranium  is  present  in  the  form  of  uranium  dioxide  that  is  20%  enriched 

235 
in  U  "  .   The  sampling  procedure  is  outlined  in  figure  1  above.   The 

results  are  described  in  the  sections  following  this  and  a  description 

of  equipment  used  is  given  in  Appendix  1. 
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Ill .   The  Gamma  Spectrum  of  Reactor  Fuel  Gas 
A.   A  Qualitative  Analysis 

The  gamma  spectra  obtained  from  the  fuel  gas  have  an  extremely 
complex  and  dynamic  nature.   The  total  spectrum  contains  about  thirty- 
four  peaks  which  can  be  shown  by  detailed  analysis  to  represent  well  over 
forty  gamma  rays.   The  previously  reported  spectra  (3)  contained  thirteen 
photopeaks.   Illustrative  fuel  gas  spectra  are  shown  in  Figures  2,  3,  4, 
and  5.   These  spectra  are  representative  of  the  approximately  400  spectra 
taken  during  this  project.   The  number  given  above  a  peak  in  a  spectrum 
corresponds  to  a  peak  number  in  Table  I,  which  contains  a  list  of  the 
corresponding  energies.   The  decade  markers  on  the  activity  ordinate  are 
useful  only  for  guaging  the  relative  heights  of  peaks  in  the  same  spectrum. 
The  natural  background  radiation  was  significant  in  Spectra  L,  M,  and  N, 
so  80%  of  the  background  was  subtracted  prior  to  graphing. 

Nearly  all  of  the  peaks  indicate  by  their  shape,  shift  in  peak 
energy,  and  complex  decay,  that  they  are  a  composite  of  several  gamma  rays. 
Obvious  examples  of  peaks  that  show  the  characteristics  mentioned  are:  the 
peak  at  0.255  Mev  which  shows  two  extremely  different  half  lives;  the 
peak(s)  that  are  initially  about  0.43  Mev  which  is  apparently  at  least 
three  peaks  of  widely  different  half  lives;  and  the  very  wide  flat-topped 
peak  with  a  nominal  energy  of  0.87  must  include  several  gamma  rays  of 
slightly  different  energies  and  very  similar  half  lives.   Only  one  marked 
qualitative  difference  was  seen  in  the  spectra  obtained  from  21  irradia- 
tions of  the  fuel  and  that  was  in  the  0.255  Mev  peak.   The  height  of  the 
peak  was  much  larger  in  the  spectra  of  one  of  the  runs.   There  was  no  readily 
apparent  reason  for  this  phenomenon,  and  the  matter  was  not  pursued.   The 
peak  above  3  Mev  was  seen  only  in  one  spectrum,  but  it  was  one  of  the  earli- 
est 4  Mev  spectra  taken  during  the  project.   The  length  of  irradiation  and 

neutron  flux  were  varied,  and  they  only  affected  the  amount  of  gas  produced. 
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Photopeaks  Present  in  Reactor  Fuel  Gas  Gamma  Spectra 


Previously 

Energy 

(Mev) 

Reported 

1 

0.012 

+0 . 004 

2 

0.032 

+0 . 002 

(3) 

3 

0.051 

+0 . 002 

4 

0.079 

+0.002 

(2),  (3) 

5 

0.148 

+0.005 

(3) 

6 

0.161 

+0.005 

(3) 

7 

0.201 

+0 . 005 

8 

0.233 

+0 . 005 

9a 

0.253 

+0.005 

(2),  (3) 

9b 

0.265 

+0.005 

10 

0.284 

+0 . 005 

11 

0.310 

+0.005 

(3) 

12 

0.364 

+0 . 005 

13a 

0.410 

+0 . 005 

(3) 

13b 

0.418 

+0 . 005 

14 

0.439 

+0.005 

15 

0.468 

+0.005 

16 

0.527 

+0.005 

17 

0.548 

+0.005 

Table  1 


Photopeaks  Present  in  Reactor  ^uel  Gas  Gamma  Spectra 


Energy 

(Mev) 

Previously 
Reported 

18 

0.608 

+0.005 

(2) 

19 

0.638 

+0.005 

20a 

0.81 

-0.92 

20b 

0.83 

-0.90 

. 

21 

0.886 

+0.005 

(3) 

22a 

1.01 

+0.01 

22b 

1.03 

+0.01 

23a 

1.14 

+0.01 

23b 

'1.17 

+0.01 

24 

1.23 

+0.01 

(3) 

25 

1.44 

+0.01 

(3) 

26 

< 

1.58 

+0.01 

(3) 

27 

1.74 

+0.01 

28 

1.85 

+0.01 

(3) 

29 

2.00 

+0.01 

30 

2.05 

+0.01 

31 

2.22 

+0.01 

32a 

2.37 

+0.01 

32b 

2.40 

+0.01 

(3) 

33 

2.59 

+0.01 

34 

3.3 

+0.2 

Table  I  (continued) 


17 


O 


5  7      14 

I  !      !  ,7 


Ulil' 


ll  ' 


22  b 


24 


9d 


X^ 


25 


27 


29 


\A.  A 


31 


,)' 


V  ^  A 


22a 

1 23a  V"ii, 

21  Stf      V 


34 


v„ 


w 


i  >  *  • 


V     t* 


32a      H^. 


V 


33 


V', 


^   iy 


20b  *\   J    t  D 

23b'        W  W    A  I 

.  26  '\A(f>  *4 

w      n  30       Vii     Si, 

^W   ■* I «I       ,     I       »ff?,C>     !  ! 

I,        i  ■  i  ii  ■ 


W,* 


1 


w   vm  fl       \       i  V  Wj 


v  ' 


11    r 


Vt  ! 


i 


/'  >i  i 


'  '  uut  t  '    I 


i       k' 


<  >f  mi  niiiimiiiiiminmiirtiiiiiKiiiiim 


M 


EV 


SPECTRUM 
A 
B 
C 
D 


RUN 

COOLING  TIME 

16 

6.40  minutes 

16 

72.50  minutes 

16 

146.60  minutes 

16 

248.00  minutes 

Early  4  Mev  Spectra 

Figure  2 

COUNTING  TIME 

1  minute 

4  minutes 

10  minutes 

20  minutes 


o 


/. 


5    7      9b 
I 

A 


14 


i''^U»ji^0ifiMii^^^^ 


,  I 


i    '  l'UlU.«MlJU>liiLmU^       ^l 


A/v  \      r"\ 


17 

I 


13d  15 


"^^^■JiliigiiUiliiul'i'l^^iiiiiJl^ 


A     ft 

•  \  v  ! 


\     Aa  g 

20a     ,H 


1      i       •  '<  .   '     ^     II 


A 


\A,/  \ 


Mev 


SPECTRUM 

RUN 

COOLING  TIME 

COUNTING  TIME 

E 

15 

11.55  minutes 

1  minute 

F 

15 

41.50  minutes 

4  minutes 

G 

15 

117.08  minutes 

4  minutes 

H 

15 

Early 

282.75  minutes 

1  Mev  Spectra 
Figure  3 

20  minutes 

19 


2     4 


o 


9a 

I 

?      A 


16 


n 


f,      i  !    7  /     \  II 


J,'     A      '    j 

V 


^^^y^^^ 


w 

%^/\j 


'jW 


Mev 


SPECTRUM 
I 
J 
K 


RUN 
21 
21 
21 


COOLING  TIME 
25.36  hours 
41.16  hours 
73.88  hours 


COUNTING  TIME 
40  minutes 
40  minutes 
100  minutes 


Intermediate  1  Mev  Spectra 
Figure  4 


20 


o 
o 


2     4 


ft :.' 


t'  ill  A 

■frl/V'     f     ?    10 


16  19 

I 


Mev 

— 1 
1 

SPECTRUM 

RUN 

COOLING  TIME 

COUNTING  TIME 

L 

21 

7.149  days 

1000  minutes 

M 

21 

12.033  days 

1000  minutes 

N 

21 

21.700  days 

Late  1  Mev  Spectra 
Figure  5 

1000  minutes 

21 


B.   The  Nuclides  Probably  Present 

The  obvious  starting  point  for  determining  what  nuclides  are 
present  in  wuch  a  complex  mixture  is  to  determine  what  nuclides  may  be 
present.   The  first  step  was  a  careful  study  of  a  standard  triangular 
chart  of  the  nuclides  and  the  result  is  shown  in  Figure  6.   The  areas 
delineated  in  Figure  6  are  all  the  decay  chains  that  may  be  present,  con- 
sistent with  the  constraints  implicit  in  the  scope  of  the  project,  and 
experimental  technique.   The  two  primary  constraints  were,  of  course, 
volatility  and  half  life.   About  1  minute  elapsed  between  ceasation  of 
irradiation  and  drawing  the  gas  into  the  sample  ampul,  and  another  3  to 
7  minutes  elapsed  before  counting  started.   This  effectively  eliminated 
the  volatile  isobaric  chains  with  atomic  masses  of:   92,  93,  94,  95,  141, 
142,  143,  and  144.   In  general,  nuclides  that  occupy  a  position  in  the 
chart  immediately  adjacent  to  the  areas  shown  terminate  the  chains  because 
they  are:   stable;  daughters  to  which  all  non-negligible  decay  is  by 
beta  (-)  to  the  ground  state;  not  a  fission  product;  or,  because  no  nuclide 
exists  at  that  position. 

The  second  step  was  to  examine  the  selected  groups  of  nuclides 
in  order  to  estimate  the  relative  abundances  of  each  decay  chain.   There 
are  three  pertinent  indicators  for  this.   The  first  is  the  familiar  plot 
of  the  Fission  Yield  curve  for  thermal  neutrons  as  a  function  of  A,  such 
as  the  one  in  Figure  7,  which  indicates  the  relative  production  of  each 
isobaric  decay  chain.   This  can  further  be  refined,  by  noting  the  actual 
Fission  Yields  that  have  been  measured  experimentally,  when  they  are 
available.   An  excellent  compilation  is  the  one  by  Katcoff  (7).   The  third 
useful  indicator  is  Independent  Yield  (8).   This  is  an  empirical  parameter 
of  each  decay  chain  that  gives  an  estimate  of  the  relative  yield  of  each 
member  of  the  decay  chain,  and  has  been  determined  for  atomic  numbers  in 
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ThermaJ  Neutron  Fission  Yield  (%)  of  U     (7) 

Figure  7 


the  ranges  of  78  to  95  and  128  to  145.   The  Independent  Yield  is  a  function 
(a  Gaussean  curve)  of  the  quantity  Z  -  Z  .   A  table  of  values  of  Z   for  appli- 
cable values  of  atomic  mass  is  given  in  Table  II.   So  for  a  nuclide  of  given 
atomic  mass  and  atomic  number,  one  determines  Z  -  Z  ,  and  the  IY  is  normally 
read  from  the  curve.   This  can  be  a  tedious  process  if  many  nuclides  are  being 
considered,  so  I  have  translated  the  graph  into  the  more  easily  readable  tabu- 
lar form,  and  this  is  shown  in  Table  III.   Finally,  I  have  summarized  all  the 
information  from  these  three  sources  in  Figure  8.   The  (m)  in  Figure  8  indi- 
cates the  number  applies  to  the  metastable  isomer.   This  information  coupled  • 
with  half  life  values  shows  that  the  following  nuclides  would  be  expected  to 
constitute  the  most  significant  portion  of  the  activity  in  the  sample:  Kr,  A: 
88,  89,  90;  Br,  A:  87,  88,  89;  Xe,  A:  135,  137,  138,  139,  140;  I,  A:  134,  135, 
136,  137,  138.   The  presence  of  most  of  these  nuclides  was  seen  in  the  gamma 

spectra  of  the  fuel  gas. 

24 


Z   vs.  A  Values  for  the  Thermal  Fission  of  U 


A 

Z 

P 

78 

31.0 

79 

31.3 

80 

31.6 

81 

32.0 

82 

32.3 

83 

32.7 

84 

33.1 

85 

33.5 

86 

33.9 

87 

34.4 

88 

34.9 

89 

35.4 

90 

35.8 

91 

36.3 

92 

36.8 

93 

37.4 

94 

37.8 

95 

38.4 

96 

38.2 

97 

38.65 

98 

39.1 

A 

Z 

P 

128 

50.2 

129 

50.3 

130 

50.4 

131 

50.8 

132 

51.3 

133 

51.5 

134 

51,8 

135 

52.4 

136 

52.6 

137 

53.3 

138 

53.4 

139 

53.8 

140 

54.3 

14X 

55.0 

142 

55.4 

143 

55.9 

144 

56.4 

145 

56.8 

146 

57.2 

147 

57.7 

148 

58.1 

149 

58.4 

150 

58.7 

Table  II 
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Z  -  Z  vs.  Independent  Yield  (%) 
P 


Independent 

-  z 

Yield 

p 

0.0 

58 

0.1 

56 

0.2 

54 

0.3 

52 

0.4 

50 

0.5 

44 

0.6 

40 

0.7 

36 

0.8 

31 

0.9 

26 

1.0 

19 

1.1 

15 

1.2 

11.5 

1.3 

10 

1.4 

8 

Table  III 


Independent 

-  z 
p 

Yield 

1.5 

6 

1.6 

3.8 

1.7 

3.0 

1.8 

2.2 

1.9 

1.5 

2.0 

0.9 

2.1 

0.6 

2.2 

0.4 

2.3 

0.25 

2.4 

0.1 

2.5 

0.08 

2.6 

0.04 

2.7 

0.02 

3.2 

0.005 
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IV.   Photopeak  Analysis 
A.   General 

The  gamma  spectra  obtained  from  the  fuel  gas  were  calibrated, 

in  the  standard  manner,  by  superimposing  on  them  the  spectra  of  several 

137 
useful  standards.   Those  most  frequently  used  were:  Cs    :  32  Kev  and 

0.662  Mev;  Au198:  69  Kev  and  0.411  Mev;  Mn   :  0.842:  and  Co      :  1.17, 

1.33,  and  2.50  Mev.   The  multichannel  analyzer  used  did  not  seem  to 

exhibit  the  linearity  of  energy  per  channel  that  it  has  been  credited 

with  (3),  however,  these  nuclides  provided  sufficient  points  to  determine 

the  photopeaks  in  the  fuel  gas  spectra  to  within  a  few  Kev,  with  good 

assurance.   The  energy  of  each  photopeak  was  calculated  by  measuring  the 

distance  along  the  energy  scale  from  the  nearest  standard  peak,  using  the 

energy  per  channel  ratio  of  the  nearest  bracketing  peaks. 

Parenthetically,  it  should  be  noted  that  the  standards  used  did 
not  determine  the  peaks  in  the  fuel  gas  spectra  above  1  Mev  with  nearly 
as  great  a  precision  as  those  below  1  Mev.   This  was  considered  to  be  of 
marginal  significance  since  the  peaks  in  this  region  were  not  prominent 
enough  to  conduct  a  detailed  half  life  analysis  using  the  methods  described 
below.   As  a  rule,  any  of  the  nuclides  of  interest  with  gammas  in  the 
higher  range  also  had  peaks  in  the  lower  range,  so  there  was  sufficient 
information  contained  in  the  latter  region  of  the  spectra. 

Now  then,  the  energy  calibration,  coupled  with  the  earlier  analysis 
of  the  possible  and  probable  nuclides,  allows  a  correlation  between  the 
nuclides  and  the  photopeaks.   One  other  tool  of  analysis  is  available, 
namely  half  life  determination. 
B.   Half  Life  Analysis 

Several  efforts  were  made  to  determine  the  short  and  intermediate 
component  half  lives  of  photopeaks  in  the  following  intervals:  0.14  -  1.17 
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Mev;  0.17  -  0.24  Mev;  0.24  -  0.29  Mev;  0.33  -  0.40  Mev;  0.40  -  0.48  Mev 
(the  members  of  the  trio  were  analyzed  both  collectively  and  individually); 
0.49  -  0.57  Mev:  and  the  photoplateau,  0.81  -  0.91  Mev.   Referring  to 
Figure  2  and  Figure  4,  these  intervals  correspond,  respectively,  to  peaks: 
5;  7;  9a  and  9b;  12:  13a,  13b,  14,  15:  16;  20a.  The  half  life  determina- 
tions were  made  by  following  the  decay  of  the  photopeaks  with  a  single 
channel  analyzer  which  gave  repetitive  outputs  of  the  integrated  number 
of  disintegrations  under  the  photopeak  per  unit  time.  These  data  were 
analyzed  using  the  computer  program,  FRANTIC  II  (9),  which  has  been 
adapted  for  use  on  the  NPS  IBM  360  Duplex  computer. 

FRANTIC  is  a  very  capable  and  flexible  program  designed  for 
analysis  of  radioactive  decay  problems  and  is  widely  used  by  chemists  and 
physicists  at  this  institution.  There  are  several  modes  in  which  it  can  be 
used,  and  it  was  used  in  a  mode  in  which  data  of  the  type  lust  described 
are  inputed,  along  with  estimates  of  the  component  half  lives.  The  pro- 
gram then  uses  a  iterative  least  squares  method  to  determine  which  of  the 
half  life  estimates  best  fit  the  input  data.  The  program  can  cope  with 
photopeaks  containing  up  to  10  components,  however,  to  get  acceptable 
results  each  component  half  life  should  differ  from  the  one  of  nearest 
magnitude  by  about  a  factor  of  two. 

FRANTIC  was  used  to  conduct  about  forty  analyses  of  the  energy 
intervals  listed  above.  Tyoically,  the  program  was  used  to  look  for  one 
to  six  components  during  each  analysis.  The  early  results,  with  only  a 
few  exceptions,  did  not  show  the  expected  number  of  comnonents.  Usuallv, 
only  one  half  could  be  found  that  agreed  well  with  accepted  literature 
values  for  the  probable  components.   It  was  believed  that  results  of  this 
type  were  obtained  for  two  reasons.  The  first  is  that  the  peaks  of  interest 
were  too  strongly  influenced  by  adjacent  peaks;  this  can  be  seen  easily  in 
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the  spectra.   A  second  complicating  effect  is  that  the  single  channel 
analyzer  measures  the  portion  of  the  Compton  scatter  distribution  under 
the  peak  of  interest  for  all  peaks  of  higher  energy.   The  effect  was  that 
these  interferences  masked  the  contributions  of  the  minor  components  in 
the  peak.   So  a  technique  was  developed  to  remove  these  two  interferences. 
C.   The  Two  Analyzer  Technique 
(1)   Description 

The  principle  is  to  follow  the  activity  decay  of  a  small  energy 
range  on  each  side  of  the  peak  being  analyzed.   Then  assume  this  general 
level  of  activity  is  constant  between  the  two  intervals,  i.e.,  under  the 
peak  of  interest,  and  subtract  the  appropriate  proportionate  amount  from  the 
activity  of  the  photopeak.   This  is  accomplished  with  two  single  channel 
analyzers  in  the  following  manner.   After  appropriate  calibration,  one 
analyzer  is  set  with  an  energy  window  that  includes  the  photopeak.   The 
second  analyzer  is  set  on  a  wider  energy  window  and  slaved  to  the  first 
analyzer  so  that  it  counts  simultaneously  with  it.   After  the  data  are  taken 
the  counts  per  counting  interval  from  the  first  scaler  are  subtracted,  point 
for  point,  from  the  same  data  from  the  second  scaler.   The  result  is  then 
proportioned  to  the  energy  interval  of  the  first  analyzer  and  subtracted, 
point  for  point,  from  the  counts  registered  by  the  first  analyzer. 

The  data  manipulation  can  be  expressed  more  simply  algebraically. 
If:  W,  =  width  of  energy  window  for  analyzer  monitoring  the  peak;  A^  = 
number  of  counts  registered  by  first  analyzer,  and  W2  and  A2  are  the  corre- 
sponding quantities  for  the  second  analyzer;  then  the  corrected  FRANTIC  input 

W-, 
data  point,  P,  equals:  A-,  -  _ — £_  (A?-A-, ) .   This  is  readily  simplified  to: 

W2~W^ 

P  =  C]Ai  -  C2A2,  where  C2  =    l        and  Ci  =  C2  +  1.   Then  it  requires  only 

W2-W! 

a  simple  modification  to  the  computer  program  to  input  A-^ ,  A2 ,  Ci  and  C2 ,  and 

have  P  computed  for  each  data  point  and  used  in  the  calculations  of  the 

program. 
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Experience  showed  that  it  was  necessary  to  determine  C.  and  C 
in  a  way  that  obviated  the  uncertainty  in  W  and  W„,  rather  than  calcu- 
lating them  directly  from  W-  and  W„.   The  constants  can  be  determined 
experimentally  by  using  a  standard  sample  that  has  a  constant  level  of 

Compton  distribution  across  the  ranges  to  be  analyzed  with  the  single 

60 
channel  analyzers.   Co   is  ideal  for  windows  below  1  Mev.   If  the  analvers 

are  set  with  the  desired  windows,  and  C  and  C~  are  correct,  then  when 

looking  at  a  level  spectrum,  the  corrected  data  point,  P,  must  equal  zero. 

Al 

So  from  the  basic  equation  for  P,  we  have  C~  =  - — _■  ■. ■  .   So  to  determine 


C  and  C  it  is  only  necessary  to  make  a  few  counts  and  calculate  C_.   One 
important  advantage  with  this  method  of  determining  C.  and  C~  is  that  their 
precision  is  limited  only  by  the  length  and  number  of  counts  one  is  willing 
to  make  to  determine  A.     and  A„. 

The  Two  Analyzer  Technique  of  data  treatment  was  tested  with  two 
test  runs  and  the  results  obtained  from  FRANTIC  are  listed  in  Table  IV 

below.   The  test  was  to  analyze  the  0.85  Mev  peak  in  a  sample  that  con- 

9  7  S6 

tained:  Mg   ,  t  1/2  =  9.1  minutes,  photon  energy  =  0.835  Mev;  Mn  ,  t  1/2  = 

64 
2.58  hours,  photon  energv  =  0.845  Mev;  and  5.1  minute  Cu   and  12.9  hr 

66 
Cu   whose  principle  spectral  peaks  are  the  0.51  Mev  and  1.02  Mev  annihi- 
lation peaks.   The  Cu  isotopes  were  added  to  provide  a  high  level  of 

27 
Compton  scatter  distribution.   The  Mg   did  not  activate  sufficientlv  to 

be  detected,  but  the  Two  Analyzer  Techniaue  succeeded  very  well  in  remov- 
ing the  contribution  made  to  the  peak  by  the  Compton  distribution  of  the 
higher  energy  peaks.   The  sample  for  the  second  run  was  intended  to  differ 
from  the  first  onlv  in  that  the  amount  of  Mn  was  about  one  one-half  of 
that  in  the  sample  used  in  the  first  run.   However,  the  magnesium  salt 
used  in  the  sample  of  the  second  run  contained  an  impuritv  that  activated 
so  well  that  it  masked  even  the  Cu  activity.   The  unidentified  impuritv 
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had  a  nominal  half  life  of  a  few  hours  and  contributed  the  35.2  minute 
half  life  shown  in  Table  IV. 

PHOTOPEAK  ANALYSIS 


Test  Pun  1 


Test  Run  2 


Photopeak  Data 
Only 

2.72  hr 
5 . 1  min 

2.74  hr 
35.2  min 


Two  Analyzer 
Technique 

2.66  hr 


2.60  hr 


Table  IV 

(2)   Application  to  Fuel  Oas  Spectrum 

The  Two  Analyzer  Technique  was  applied  to  the  photopeaks  listed 
at  the  beginning  of  Section  IV  R  above  and  it  gave  improved  results,  in 
several  cases.   The  results  were  improved  in  the  same  manner  as  that 
shown  in  the  test  runs,  i.e.,  it  eliminated  spurious  half  lives  that  were 
contributed  by  the  Compton  distribution.   However,  it  did  not  allow  detec- 
tion of  additional  minor  component  half  lives  in  the  photopeak,  as  had  been 
hoped  for.   For  a  comparison  of  the  results  obtained  from  applying  both 
methods  to  fuel  gas  spectra,  see  the  entries  in  Table  V  below  for  Peak  No.  7 
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V.   Results 

The  heart  of  the  experimental  results  has  been  displayed  in  the 
soectra  of  Figures  2-5.   Because  of  the  many  peaks  involved,  and  the  even 
greater  amount  of  nuclidic  data  that  must  be  correlated,  the  interpreta- 
tion of  the  spectra  is  done  largely  in  tabular  form  below.   The  peak 
numbers  listed  in  Table  V  correspond  to  the  peak  numbers  given  in  Table  I 
and  Figures  2-5.   The  information  on  the  nuclides  was  taken  largely  from 
the  very  convenient  compilation  by  Crouthamel  (10) .   A  (a)  after  a  half 
life  value  indicates  it  was  determined  using  photopeak  data  only.  A  (b) 
indicates  it  was  determined  using  the  Two  Analyzer  Technique.  A  number  in 
parenthesis  represents  a  bibliography  entry,  as  usual.  Where  several 
acceptable  measurements  were  obtained  all  are  listed.   FRANTIC  outputs  a 
standard  deviation  for  each  half  life  determined,  and  this  is  not  listed 
in  order  to  conserve  space,  and  because  the  FRANTIC  result  is  not  shown  if 
the  standard  deviation  was  more  than  5%  of  the  half  life  literature  value. 
Photon  energy  is  stated  in  Mev.   Where  a  nuclide  assigned  to  a  peak  is 
enclosed  in  parenthesis,  the  parenthesis  indicate  that  such  an  assignment 
is  possible  but  not  well  founded.   A  list  of  the  nuclides  discovered  is 
shown  in  Table  VI. 
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SUMMARY  OF  EXPERIMENTAL  RESULTS 


Peak 
No. 


Half  Lives 


1.    See  Figure  A 


Possible 
Nuclides 


v   85 
Kr 

11Y 

B   91 

Sr 

9.7h 


82 
Br 

35.  3h 

Y91m 
50m 


Nuclide 
Assignment /Remarks 

Peak  only  displayed  in  Spectrum  J, 
but  it  occurred  in  late  spectra  at 
least  6  different  irradiations.  It 
was  deleted  from  Spectra  H  through 
L  during  printing. 

(Kr85>  (Br82)  (Sr91=  Y91) 


2. 


21h 

growth  (3) 


All  of  Br  and  Kr  No  growth  seen  in  any  spectra  of 

nuclides  and     this  project. 

their  daughter 

products. 

See  Figure  6. 


3. 


5,  6. 


133 
Xe 

0.081, 

0.053 

5.28d 

5.29d  (3) 

Xe133 

131 

0.080 
8.09d 

A.69h, 
21. Ah  (a) 

5.  Kr85m 
0.150 
4.36h 

v   88 
Kr 

0.166 

2.77h 

3.9Ah  (a) 

Xe135 

A.97h  (b) 

0.1A0 
0.2h 

2 .  3h , 

7.9Ah, 

3.6m, 

21  .9m  (a) 

6.  Xe"1 
0.16A 
1  2d 

Xe 


133 


Xe 


133 


Kr85m,  Kr88,  (Xe135) 


Xe 


131 


4.5h  (3) 


Table  \' 
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SUMMARY  OF  EXPERIMENTAL  RESULTS 


Peak 

Poss: 

Lble 

No. 

Half  Lives 

Nucl: 

Ides 

7. 

25.9m, 

v   88 
Kr 

*  89 
Kr 

1.2h, 

0.191 

.  0.22 

3.1h  (a) 

2.77h 

3.2m 

2.51h  (b) 

136 

0.210 

3. Oh  (3) 

1.4m 

8. 

See  Figure  5 

v   133m 
Xe 

0.233 

2.3d 

b. 

1.7m,  30. 5m, 

Xe135 

i136 

11. 5h  (a) 

0.250 

0.270 

9.13 

1.4m 

11.03h  (a) 
10.75h  (a) 
10.65h  (a) 
8.58h  (b) 
9.1h  (3) 
2.5d  (3) 


Kr85m 


See  Figure  5 

i131 

i136 

0.284 

0.270 

8.09d 

1.4m 

See  Figure  4 

v   85m 
Kr 

«  140 
Ba 

4.2h  (3) 

0.305 

0.301 

4.36h 

12. 8d 

See  Figure  5 

Xe135 

.131 

0.36 

0.364 

9.13 

8.09d 

Table  V 

(continued) 
35 

Nuclide 
Assignment /Remarks 


v   88 
Kr 


.  133m 
Xe 


_  135   ,_136, 
Xe   ,  (I   ) 


I131 


I131 


SUMMARY  OF  EXPERIMENTAL  RESULTS 


Peak 

Poss 

ible 

Nuclide 

No. 

Half  Lives 

Nucl 

ides 

Assignments /Remarks 

13a, 

13,  14,  15: 

I136 

Kr87 

_  87   _  138    138 
Kr   ,  Cs    ,  Xe 

13b, 
14, 

4.6m,  27.7m, 
1.56h  (a) 

0.388 
1.4m 

0.403 
1.3h 

(Br8Am)  (Xe137) 

15. 

13:34m,  1.7h, 

Xe138 

.  140 
La 

6.5h  (a) 

0.42 
17m 

0.438 
40.  2h 

13:4. 5h 

2.2h  (3) 

Br 

0.44 

6.0m 

Xe137 

0.44 

4.2m 

Cs138 

0.463 
32.2m 

(5) 

16 


See  Figure  5 


I133 

0.53 
20. 8h 

I131 
0.51 
8.09d 


..  135m 
Xe 

0.53 

15.6m 


.131 


(I133) 


17.    See  Figure  2 


Sr91 

9.7h 


Cs138 
0.55 
32  min 


y91m 

0.55 
50.3m 

i,  82 
Br 

0.554 

36  h 


Cs 


138 


18 


See  Figure  4 


Xe135 

0.61 

9.13h 


Xe 


135 


19.    See  Figure  5 


i131 

0 .  637 

8.09d 


Sr91 
0.645 

(>.  /h 


.131 


Table  V  (continued) 
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SUMMARY  OF  EXPERIMENTAL  RFSULTS 


Peak 

Possible 

Nuclide 

No. 

Half  Lives 

Nuclides 

Assignments /Remarks 

20a, 

40.7m, 

x132    I133 

20b, 

2.91h  (a) 

0.78    0.85 

21. 

2.3h    21h 

2.54h, 
2.68h  (b) 

T134    „  87 
I       Kr 

0.85    0.85 

2.7h  (3) 

52.5m   76m 

v   88    „  84m 
Kr      Br 

0.85    0.88 

2.77h   6.0m 

y89m    Rb88 

0.915   0.908 
50d     17.8m 

B  84    _  91 
Br      Sr 

0.81    0.93 

32m     9.7h 

22. 

About  20  minutes, 

Kr89    RbM 

(Kr89-Rb89)  and 

initially. 

3.2m    1.05 

15.4m 

T135    _  84 
I       Br 

1.04    1.01 

6.68h   32m 

/T135    _  84N 
(I    or  Br  ) 

23a, 

134     135 

,134     135, 
(I    or  I   ) 

23b. 

1.10    1.14 
52.5m   6.68h 

24. 

2.7h  (3) 

rr   89    nl_89 
Kr      Rb 

3.2m    1.26 

15.4m 

(Kr89-  Rb89) 

Table  V  (continued) 
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SUMMARY  OF  EXPERIMENTAL  RESULTS 


Peak 
No. 


25. 


Half  Lives 


Possible 

Nucl 

ides 

132 

133 

1.4 

1.4 

2.26h 

21h 

Cs138 

Sr91 

1.43 

1.4 

32.2m 

9.7h 

R  84 
Br 

1.46 

6m 

Nuclides 
Assignments /Remarks 


26.    2.7h  (3) 


v   88 
Kr 

1.55 

2.8h 

.  140 
La 

1.60 

40.  2h 


Rb89 
1.55 
15m 


Kr 


88 


27,  28.   2.7h  (3) 


I135 

v   88 
Kr 

v   88 
Kr   : 

Rb88 

1.72 

2.87h 

6.68h 

Rb88 

v   87 
Kr 

1.85 

1.75 

17.8m 

1.3h 

29.    20  minutes; 
See  Figure  2 


138 


Xe 

2.01 

L7m 


Xe 


138 


30, 


v   87 
Kr 

2.05 

1.3h 


Kr 


87 


Table  V  (continued) 
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SUMMARY  OF  EXPERIMENTAL  RESULTS 


Peak 

Half 

Lives 

Poss 

ible 

Nuclides 

No. 

Nucl 

ides 

Assignments /Remarks 

31 

Xe138 

2.2 

17m 

Kr 
3.2m 

„  138 
Cs 

2.21 

32.2m 

Rb89 
2.20 
15.4m 

Xe138,  Cs138 

32a, 

2.9h 

(3) 

_  88 
Kr 

Br 

Kr88  (Br84) 

32b. 

2.40 
2.77h 

2.47 
32m 

33. 

V   87 
Kr 

2.57 

1.3h 

Cs138 

2.63 

32.2m 

™_89 
Rb 

2.59 

15.4m 

(Rb89)   (Cs138) 

34. 

Very  s 

hort. 

Rb89 

Br 

See  Figure  1. 

3.0, 

3.03, 

3.5 

3.91 

15m 

32min 

Table  V  (continued) 
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NUCLIDFS  IDFNTIFIED 


Volatile  Nuclides: 
,131 


v   131  v  133m  v   133  v  135  v  138 

Xe    ,  Xe    ,  Xe    ,  Xe    ,  Xe 

_  85m  _  81    88 

Kr    ,  Kr   ,  Kr 


Active  Daughters: 

Cs138 
Rb88 

Nuclides  Probably  Present: 


i133  (ii),  i134,  i135  (id,  i136 

Xe137 

Kr88,  Kr89 
Br84"1,  Br84 
Rb89 


Table  VI 
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VI .   Conclusion 

The  dramatic  increase  in  the  complexity  of  the  fuel  gas  spectra 
over  earlier  core  gas  spectra,  and  the  concomitant  increase  in  the  amount 
of  nuclidic  information  contained  therein,  certainly  vindicated  the  new 
method  of  obtaining  the  gas. 

The  more  complex  spectra,  however,  require  new  and  more  capable 
automated  and  computerized  means  of  half  life  determination,  if  all  of 
the  information  available  in  the  new  spectra  is  to  be  retrieved.   One 
approach  that  should  be  considered  is  to  conduct  a  detailed  mathematical 
analysis  of  the  various  equilibria  occurring  between  active  parents  and 
active  daughters.   Alternatively,  the  gases  could  be  separated  rapidly 
prior  to  applying  gamma-ray  spectrometry. 

The  Two  Analyzer  Technique  of  half  life  determination  is  a  valuable 
improvement  in  the  former  direction,  but  it  apparently  was  not  sufficiently 
capable  to  cope  with  the  fuel  gas  spectrum.   The  latter  route  has  been 
tried  (11)  in  work  carried  on  concurrently  with  this  project,  and  the 
fuel  gas  was  fractionated  chromatographically  into  iodine,  krypton,  and 
xenon  fractions.   The  resulting  spectra  were  much  less  complex,  but  the 
process,  in  its  current  state  of  development,  is  so  slow  that  short  lived 
nuclides  cannot  be  detected. 
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Appendix  I 


EQUIPMENT 


1.  Counting  Chamber:   A  2-inch  thick  lead  cube,  22  inches  square  by 

28  inches  high,  lined  with  1/8-inch  aluminum 
and  1/4-inch  high  copper  content  brass. 

2.  Scintillation  Detector:   Harshaw  12SW12-W3,  3-inch  by  3-inch  cylin- 

drical well  type  Mai  (Th)  crystal  with  matched 
photomultiplier  tube. 

3.  512  channel  multi-channel  analyzer,  consisting  of  the  following 
components: 


(a 
(b 
(c 
(d 
(e 
(f 
(g 


Nuclear  Data,  Inc  ND-180E   Analog  to  Digital  Converter 

Nuclear  Data,  Inc  ND-180M   512  Channel  Memory  Unit 

Nuclear  Data,  Inc  ND-180R   Read  out  Control  Unit 

Textronix,  Inc,  Rm  503  Oscilloscope 

Tally  Corp,  Model  1506  Reader/Perforator 

Teletype  Corp,   Teletypewriter 

Hewlitt-Packard,  Co.,  Mosley  7590  CR  Automatic  Plotting  System 


4.  Single  channel  analyzer  (2  ea) :   Hewlitt-Packard  Model  5201L  Sealer- 
Timer 

5.  Rapid  digital  printer  (2  ea) :   Hewlitt-Packard,  Co.,  562  A  Digital 
Recorder. 
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